HU-EP-09/03 



Pion condensation in the Gross-Neveu model with nonzero baryon and isospin 

chemical potentials 

D. Ebert ^ and K. G. Klimenko ^ 

^ Institute of Physics, Humboldt- University Berlin, 10115 Berlin, Germany 
^ Institute for High Energy Physics, 14-2281 Protvino, Moscow Region, Russia 
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During the last decade great attention was paid to the investigation of the QCD phase diagram in terms of baryonic 
as well as isotopic (isospin) chemical potentials. First of all, this interest is motivated by experiments on heavy-ion 
collisions, where we have to deal with dense baryonic matter which has an evident isospin asymmetry, i.e. different 
^. neutron and proton contents of initial ions. Moreover, the dense hadronic/quark matter inside compact stars is 
I ' also isotopically asymmetric. Generally speaking, one of the important QCD applications is just to describe the 
dense and hot baryonic matter. However, in the above mentioned realistic situations the density is rather small, and 
weak coupling QCD analysis is not applicable. So, different nonperturbative methods or effective theories such as 
chiral effective Lagrangians and especially Nambu - Jona-Lasinio (NJL) type models [l| are usually employed for the 
consideration of the properties of dense and hot baryonic matter under heavy-ion experimental and/or compact star 
conditions, i.e. in the presence of such external conditions as t emp erature and chemical potentials, magnetic field, 
finite size effects etc (see, e.g., the papers 0, H, |3, H, H, 0, IM Hi [l0| and references therein). In particular, the color 
superconductivity fl, [H| as well as charged pion condensation [ill. [l2l. fT3. fl3. [Tsl. [iGj phenomena of dense quark matter 
QQ , were investigated in the framework of these QCD-like effective models. 
T— I ' It is necessary to note that an effective description of QCD in terms of NJL models, i.e. through an employment 
, of four-fermionic theories in (3-|-l)-dimensional spacetime, is usually valid only at rather low energies and densities. 
Besides, at present time there is the consensus that another class of theories, the set of (l+l)-dimensional Gross- 
^\ Neveu (GN) type models [i3i[3j '^^^ ^^^^ be used for a good qualitative consideration of the QCD properties without 
any restrictions on the energy/density values, which is in an encouranging contrast with NJL models. Indeed, the 
^ GN type models are renormalizable, the asymptotic freedom and spontaneous chiral symmetry breaking are another 
properties inherent both for QCD and GN theories etc. In addition, the fJ. — T phase diagram is qualitatively the same 
in QCD and GN model [l^, [lo, Hi], [13, HI] (here fi is the quark number chemical potential and T is the temperature). 
The GN type models are also very suitable for the description of physics in such quasi one-dimensional condensed 
matter systems as polyacetylene [2J] . Due to the relative simplicity of GN models in the leading order of the large 
A^c-expansion (Nc is the number of colored quarks), their use is very convenient for the application of nonperturbative 
methods in quantum field theory [25j . Moreover, it is necessary to note that just in the leading order of the large 
iVc-expansion the well known no-go theorem by Mermin-Wagner-Coleman [26j apparently forbiddin g th e spontaneous 
breaking of continuous symmetries in the considered (H-l)-dimensional models is not valid [2ll. [2^. l23j . 

Thus, such phenomena of dense QCD as color superconductivity, where the color group is broken spontaneously, 
and charged pion condensation, where the spontaneous breaking of the continuous isospin symmetry takes place. 
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might be simulated in terms of simpler (l+l)-dimensional GN type models (see, e.g., [22[ and |27[, correspondingly). 

In our previous paper ^2T\ the phase diagram of the (l-l-l)-dimensional GN model with two massless quark flavors 
was investigated under the constraint that quark matter occupies a finite space volume (see also the relevant papers 
[1^). In particular, there we have studied in the large iVc-limit the charged pion condensation phenomenon in cold 
quark matter with zero baryonic density, i.e. at = 0, but nonzero isotopic density, i.e. with nonzero isospin 
chemical potential /x/. In contrast, in the present paper in the leading order of the 1/iVc-expansion we consider the 
phase portrait of the above mentioned GN model in a more general case, where both isospin- and quark number 
chemical potentials are nonzero, i.e. /i/ 7^ and /i ^ 0. Moreover, we deal with nonzero temperature and with 
spacetime of the usual topology, x B.^. We hope that our investigations will shed some new light on the physics 
of dense and hot isotopically asymmetric quark matter which can be observed in heavy-ion collision experiments or 
might exist in compact stars, where baryon density is obviously nonzero (i.e. /x 0). Our consideration is based, 
for simplicity, on the approach with homogeneous condensates (an extension to inhomogeneous condensates was also 



2 



considered in [H, [2§|). 

The paper is organized as follows. In Section II the thermodynamic potential of the two-flavored Gross- Neveu 
model is obtained in the presence of quark number as well as isotopic chemical potentials. In Section III and IV the 
phase structure of the model is investigated at zero and nonzero temperatures, correspondingly. It turns out that at 
zero temperature, /i/ 7^ and rather small values of /i the gapped pion condensed phase occurs. However, at larger 
values of the dense quark matter phase with chirally broken phase and gapless quasiparticles is realized. 



II. THE MODEL AND ITS THERMODYNAMIC POTENTIAL 



We consider a (l+l)-dimensional model which describes dense quark matter with two massless quark flavors {u and 
d quarks). Its Lagrangian has the form 



L 



rid. 



a , f^i 
+ Y^37 



G 



{qqf + (qij^Tqf 



(1) 



where the quark field q{x) = qia{x) is a flavor doublet {i = 1,2 or i — u,d) and color iVc-plet (a = l,...,Nc) as 
well as a two-component Dirac spinor (the summation in ([T]) over flavor, color, and spinor indices is implied); 
(fc — 1, 2, 3) are Pauli matrices; the quark number chemical potential /i in ([T]) is responsible for the nonzero baryonic 
density of quark matter, whereas the isospin chemical potential fij is taken into account in order to study properties 
of quark matter at nonzero isospin densities (in this case the densities of u and d quarks are different). Evidently, 
the model ([T]) is a generalization of the (H-l)-dimensional Gross-Neveu model [l3| with a single massless quark color 
iVc-plet to the case of two quark flavors and additional chemical potentials. As a result, in the considered case we 
have a modified chiral symmetry group. Indeed, at = apart from the color SU(iVc) symmetry, the Lagrangian 
(fT|) is invariant under transformations from the chiral SUl{'2) x SUr{2) group. However, at /i/ 7^ this symmetry is 
reduced to f//3i(l) x Ui^r{1), where I3 = T3/2 is the third component of the isospin operator (here and above the 
subscripts L,R mean that the corresponding group acts only on left, right handed spinors, respectively). Evidently, 
this symmetry can also be presented as Ui^{l) x [7^/3(1), where Ui^{l) is the isospin subgroup and [7^/3(1) is the 
axial isospin subgroup. Quarks are transformed under these subgroups as q — > exp(iar3)q and q exp(ia7^r3)g, 
respectively. 

The linearized version of the Lagrangian ([1]), which contains composite bosonic fields a{x) and 7ra(a;) (a = 1,2,3), 
has the following form: 



L 



4G 



aa + TTaTTfl 



From the Lagrangian ([2]) one gets the following constraint equations for the bosonic fields 

G G 

<^{x) = -2 — (gg); na{x) = -2 — {qij^Taq). 



(2) 



(3) 



Obviously, the Lagrangian ([2]) is equivalent to the Lagrangian ([T]) when using the constraint equations Fur- 
thermore, it is clear from ([3]) that the bosonic fields transform under the isospin Ui^{l) and axial isospin UAisi^) 
subgroups in the following manner: 



?7/3(l) : a a: TT^ ^ ir^; tti — s- cos(2Q;)7ri -I- sin(2a)7r2; TT2 
{7^/3(1): T^i ^ T^i', ^ T^2] (T — > cos(2a)cr + sin(2Q;)7r3; 



cos(2Q!)7r2 — sin(2Q!)7ri, 
cos(2a)7r3 — sin(2Q!)cr. 



(4) 



Due to the transformation rules ([4]) such quantity as the thermodynamic potential (TDP) depends effectively only 
on the two combinations (cr^ -I- 7r|) and (tt^ + ""i) of the bosonic fields, which are invariants with respect to the 
[7/3(1) x [7^/3(1) group. In this case, without loss of generality, one can put tt2 — t^s ~ 0, and study the TDP as 
a function of only two variables, M = a and A = tti. Throughout the paper we suppose for simplicity that the 
condensates M and A are homogeneous quantities which do not depend on the space coordinate. In order to avoid 
the no-go theorem [2^ , which forbids the spontaneous breaking of continuous symmetries in the considered case of 
one space direction, all our considerations are performed in the leading order of the 1/iVc-expansion, i.e. in the limit 
iVc 00, where the TDP of the model looks like (see, e.g., [13]) 



AG 



^ In { [ipo + f^f [E+r] [{po + m)2 - [E^r] }, 



(5) 



where E^ = y^{E^y+A^, E^ ^ E ±v, v ^ ^//2 and E = y^pj + AP. It is clear that the TDP f7^,^(Af, A) is 
symmetric under the transformations fi — ^ — /i and/or v —> ~v. So it is sufficient to consider only the region /i > 0, 
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v > 0. Taking into account this constraint as well as integrating in ([5]) over po, we obtain for the TDP of the system 
at zero temperature the following expression: 



AG 



dpi ( 
2n I 



E+ + E^ + i^l- E+)e{ii - E+) + - 



E-^)e{^i 



e: 



;)}, 



(6) 



where 9{x) is the Heaviside step function. Since we are going to study the phase diagram of the initial GN model, 
the system of gap equations is needed: 



dM 2G 



dpi ^e{E+ 

2ttE\ 



fi)E+ , 0{E^-t,)E- 



E-, 
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dn^^,{M, A) ^ 

5A 



dp^^ e{E+-y) 0{E^-fi) ^ 

2^1 Et E^ r 



(7) 



Evidently, the coordinates Af and A of the global minimum point of the TDP ^ supply us with two order parameters 
(gaps), which arc proportional to the ground state expectation values (qq) and (qij^Tiq), respectively. If only the gap 
M is nonzero, then in the ground state of the model the axial isospin symmetry [7^/3(1) (at /i/ ^ 0) is spontaneously 
broken down. However, if only the gap A 7^ 0, then the ground state describes the phase with charged pion conden- 
sation, where the isospin Ui^{l) symmetry is spontaneously broken. Note that in this phase the space parity is also 
spontaneously broken. 



III. PHASE PORTRAIT AT ZERO TEMPERATURE 



First of all, let us consider the phase portrait of the initial GN model (1) using as a starting point the TDP ([6]). 
Since it is an ultraviolet divergent quantity, one should renormalize it, using a special dependence of the bare coupling 
constant G = G{A) on the cutoff parameter A (A restricts the integration region in the divergent integrals, \pi\ < A). 
In our previous paper [271] the following prescription for the bare coupling constant G{A) was used. 



1 



2G(A) 




(8) 



where Mq is the dynamically generated quark mass in the vacuum, i.e. at /i = and /i/ = (see below). Then, 
introducing the quantity ri^,,y(M, A; A), 



r!^^„(M,A;A) = 



AP + A2 



4G(A) 

it is possible to obtain the renormalizcd (finite) expression for the TDP: 

r!^,,(A/,A) = lim r!,,,,(Af,A;A) 

A—* 00 



/I ^{^^ + E^ + ifi- E+)9ifi - E+) + (^ - E^)e{^l -E^)] + ^, (9) 



(10) 



(The renormalizcd expression for the gap equations is obtained in the limit A — s- 00, if the replacements G G(A) 
and \pi\ < A are done in ([7|), or by a direct differentiation of the expression (fTO| .) In particular, at /i = and /// = 
we have from (flOl): 



n^^,{M, A) 



A^2 



= K)(VM2 + A2) = 



In 



+ A^ 



- 1 



(11) 



Since for a strongly interacting system the space parity is expected to be a conserved quantity in the vacuum, we put 
A equal to zero in (jlip . As a result, the global minimum of the TDP (jlip (usually, this quantity is called effective 
potential in the vacuum) lies in the point A'l — AIq, which means that in the vacuum the dynamically generated quark 
mass is just the parameter Afo introduced in ([5]). However, in the general case, i.e. at nonzero values of the chemical 
potentials, this quantity depends certainly on and obeys the gap equations ([7]). 

Numerical investigations show that local minima of the TDP (|10p can occur only on the Af- or A-axes (other 
solutions of the gap system ([7]) correspond to saddle points of the TDP PH]) ). So it is enough to restrict the 
consideration of the TDP ^ to the regions A/ = 0, A ^ (A-axis) or M 7^ 0, A (M-axis). Moreover, since the 
TDP is an even function with respect to the transformations Af — ^ — Af or A ^ — A, we will suppose further that 
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FIG. 1: The (/i, i^) phase portrait of the model at T = 
and V > Q. Here jl = v — and Mo is the quark 

mass in the vacuum. I^umber 1 cienotes the symmetric 
phase with massless quarks, number 2 - the normal quark 
matter phase with massive quarks, and 3 denotes the pion 
condensed phase. 



FIG. 2; The {t, fi) phase portrait of the model at /i/ = 
0.2Mo {t = -^,fj' = 'm^)- The solid and dotted lines are 
the critical curves of the 1st and 2nd order phase transi- 
tions, correspondingly. They are divided by the tricritical 
point A. Other notions are the same as in Fig. 1. 



M, A > 0. As a result we have 



6i(/i - A) 



A^ In 



A 



(12) 



r2^,,(M,A = 0) = Vb(Af)-H 



9{fi + u- M) 



27r 



In 



M 



Af2 



\y^~v\-M) 
2tt 



A//2 In 



\^i-u\ + ^{^i-uY-^p 



M 



\li - yy {y. - yf - NP] (l - e{M - v)B{v - y)) , 



(13) 



where the function Vq{x) is presented in pT|) . Compairing the global minima of the TDPs ((T^ - P^ . it is possible 
to obtain the global minimum point of the genuine TDP PH]) of the model. Clearly, its form (and, as a result, the 
phase of the model) depends on the values of /i and v. So, the total information about the behaviour of the global 
minimum point vs /i, can be presented by the phase portrait given in Fig. 1. (Note, it is valid only for v > 
values. In the z/ = case the phase structure is described below). There, in the corresponding (i/, /i)-plane (excluding 
the /i-axis), where p, = v = — jj^, you can see three respective phase regions denoted by the numbers 1, 



2, 3. For the values of (i',/!) from the region 1, 2, and 3 the global minimum point of the TDP pOj) has the form 
(A/ = 0, A = 0), {M ^ 0, A = 0), and {M = 0, A ^ 0), correspondingly. As a result, in the region 1 the chirally 
^7/3(1) X J7a73 (l)-symmetric phase with massless quarks is arranged. In the region 2, where the order parameter M 
is nonzero, this symmetry is spontaneously broken down to the isospin Ui^{l) subgroup (in this region the order 
parameter M is a smooth function vs /i and ly). We call this phase the normal quark matter phase, since here quarks 
dynamically acquire a mass which is equal to the order parameter M , and space parity is not broken. Finally, the 
region 3 corresponds to the charged pion condensed phase, because of the nonzero order parameter A. Furthermore, 
for all points of this region A = Mq. Note that on the phase boundaries phase transitions of the first order occur. 

It is worth to remark that in early papers (see e.g. [3, [l^]), the phase structure of the GN model was studied in the 
exceptional case of = 0, 7^ with homogeneous condensates. As it was shown there, at /i > the symmetric 

phase 1 is observed. However, at y < A/o/\/2 the normal quark matter phase is realized with a constant quark mass 
equal to A/q. ^ Let us denote this particular normal phase as phase 4, which is not depicted in Fig. 1. ^ There are 



^ If an inhomogeneous ansatz for condensates is taken into account, then in some GN phase diagrams an additional crystalline phase 
might appear at 1^ = and intermediate values of fi |23l |29|| . The existence of the crystalline phase in the general case with u ^ is not 
investigated up to now. 

We would like to stress once more that in this figure the phase 3 does not occupy the part of the /i-axis (P = 0, < A < I/V2), where 
actually the phase 4 is arranged. 
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two essential differences between phases 2 and 4. The first one manifests itself in their dynamical properties such as 
quasiparticle excitations. Recall that in the most general case the energy spectrum of the u-, d-, u-, d-quasiparticles 
(quarks) has the following form (see, e.g., [H]): 

EM^E^-f^, Edipi) = E+ - fi, E^ip,) = E+ + ^i, EM = E^+fi. (14) 

It is easily seen from ((T?| that in the phase 4, where i^ — 0, A = 0, M — Mq and fi < Afo/"\/2, there is a gap in the 
quasiparticle energy spectrum, i.e. in the phase 4 all quasiparticles are gapped excitations. In contrast, in the phase 
2 (u ^ 0) the situation is quite different and here in the energy spectrum of the M-quasiparticles the gap is absent. It 
means that for each point of the region 2 in Fig. 1 there exists a space momentum such that i?„(p*) = 0. (For 
example, the point {i> = 0.7, /2 « 0.92) is in the phase 2 (near its boundary), where A = and M « O.SlAfg. In this 
case pi « 1.59Afo etc.) So, in analogy with gapped and gapless color superconductivity ,4], the phase 2 may be called 
gapless normal quark matter phase. 

The second difference between phases 2 and 4 is a thermodynamic one and related with the density of quarks in 
both phases. It is clear that in the phase 4, i.e. at = 0, M = Mq and fi < Mo/V2, the TDP from (|13p does not 
depend on the quark chemical potential fi. As a result, we see that in this phase the quark density = —d^/dji is 
equal to zero. However, in the phase 2 the density of quarks is not zero. So the consideration of dense quark matter 
in the framework of the GN model is more adequate in terms of both nonzero isotopic and quark number chemical 
potentials, v ^ Q and n ^ 0, than in the simpler case with = 0, 7^ 0. 

It is seen from Fig. 1 that the symmetric phase 1, where M = and A = 0, is a disconnected manifold consisting of 
two different parts. (Note that in the ground state of this phase the quark number density is not zero.) The first part 
is arranged above the phase 2 and corresponds to fj, > v, whereas in the second one, which is below the phase 2, we 
have n < V. It turns out that the dynamical properties of these two regions are different. Indeed, for the part of the 
phase 1 which is above the phase 2 both u- and c?-quasiparticles are gapless. However, in the rest of the phase 1, which 
corresponds to fj, < v, only u-quasiparticles are gapless, as it is clear from the corresponding dispersion relations of 
these quasiparticles, Eu{pi) — \pi — v\ — ^j. and Ed{pi) — pi + v — ^l, in the phase 1. Hence, such dynamical properties 
of dense quark matter as transport phenomena (e.g., conductivities etc) can occur in a qualitatively different way in 
these subphases of the symmetric phase 1 . 

Finally, note that all quasiparticle excitations of the pion condensed phase 3 are gapped ones. Indeed, it is clear 
from (HH) that at A = Mq, M = and fx < Mo/y/2 the quantities (fTi|) do not vanish. So there is a gap in the energy 
spectrum of quasiparticles in the pion condensed phase 3. 



IV. PHASE PORTRAIT AT NONZERO TEMPERATURES 



Now let us stydy the influence of the temperature T on the phase structure of the considered GN model with two 
chemical potentials fi and v = /i//2. To get the corresponding thermodynamic potential fi^ 3,(Af, A), one can simply 
start from the expression for the TDP at zero temperature © and perform the following standard replacements: 



dpo 
2tt 



PO POn 



IUJ„ 



iTTT{2n+l) 



0,±1,±2, 



(15) 



n— — 00 



i.e. the po-integration should be replaced in favour of the summation over an infinite set of Matsubara frequencies 
Un- Summing in the obtained expression over Matsubara frequencies (the corresponding technique is presented, e.g., 
in [s^l), one can find for the TDP: 



fl,fJ.I,T 



{M, A) 



r '^\e+ + i?- + Tin [1 + e-''(^i-^)] + Tin [1 + e-'3(^^l+'^)] 

J -(XI ^TT L 



Tin [l + e"^(^A"^)] +Tln[l- 



(16) 



where (3 = 1/T. As in the case with zero temperature, the renormalized expression for ri^,,y_r(M, A) can be obtained 
by the replacement G G{A) (see formula (|8])) which, along with the cutting of the integration region in (|16p . 
\pi\ < A, leads in the limit A ^ 00 to the finite expression for the TDP. Numerical investigations show that all 
possible local minima of the obtained TDP lie on the M- and A-axis. So it is sufficient to deal with corresponding 
restrictions of the TDP on these axes, i.e. with the following functions. 



%,^.,AM = 0,A) 



yo(A) 

TT TT 









"1 + 











r!^,^,,^(M,A = 0) = T/o(Af)-^ 
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TT 



dpi In I 
dpi In I 



1 + e"'^(^+''"''' 



-l3(E+v+ti.) 



(17) 



(18) 
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FIG. 3: The (t, fl) phase portrait of the model at /i/ = FIG. 4: The [t, fl) phase portrait of the model at = 

1.2Mo {t = -p^, jl — jj^)- The points A, B, C are tri- 2Mo (t = -jg^, fl — jj^)- Other notions are the same as 

critical points. Other notions are the same as in Figs. 1, in the previous figures. 
2. 



where the effective potential Vo{x) is given in (jlip. E ~ ^Jj^^^^M^ ^ and E = y^pf + A^ . Comparing the global 
minima of the functions (|17p and (|18p , it is possible to establish the global minimum point of the renormalized TDP 
(fT6|) . Then, the dependence of the global minimum point vs T, /x, v defines the phase structure of the model. 

Using this prescription in our numerical investigations of the TDPs (fT7)) - p^ . we have found the three (/i, r)-phase 
portraits of the initial GN model, depicted in Figs. 2, 3, 4, for three qualitatively different values of the isospin chemical 
potentials, fij — 0.2Afo, m — 1.2Mo, and fij — 2Mo, respectively. There, the soHd and dotted fines correspond to 
curves of first- and second order phase transitions. Moreover, there are several tricritical points. A, B, C, in these 
phase diagrams. ( A point of the phase diagram is called a tricritical one, if an arbitrary small vicinity of it contains 
both first- and second order phase transition points.) 



V. SUMMARY AND CONCLUSIONS 



Recent investigations of the phase diagram of isotopically asymmetric dense quark matter in terms of NJL models 
show that their pion condensation content is not yet fully understood. Indeed, the number of the charged pion 
condensation phases of the phase diagram depends strictly on the parametrization set of the NJL model. It means 
that for different values of the coupling constant, cutoff parameter etc just the same NJL model predicts different 
numbers of pion condensation phases of quark matter both with or without electric neutrality constraint (see, e.g., 
[Hill). So to obtain a more objective information about the pion condensation phenomenon of dense quark matter 
it is important to invoke alternative approaches. One of them, which qualitatively quite successfully imitates the 
QCD properties (see also the Introduction), is based on the consideration of this phenomenon in terms of the leading 
order of the large -/Vc-technique in the framework of (l-l-l)-dimensional GN models. 

In the present paper we have studied the phase structure of the GN model (1) in terms of temperature, quark 
number (fi)- as well as isospin (ni) chemical potentials in the limit Nc — > oo. Firstly, we have found that at T = 
the charged pion condensation phase of the GN model is realized inside the chemical potential region fij > 0, and 
/i is not greater than (see Fig. 1). The corresponding quasiparticle excitations have a gap in the energy 

spectrum. In contrast, in the NJL phase diagram the pion condensation phases occupy a compact region and for some 
parameterization schemes the gapless pion condensation might occur 

Secondly, at /i/ > and rather large values of the quark number chemical potential /i we have found the normal 
quark matter phase 2 (see Fig. 1) in which chiral symmetry is broken spontaneously. Moreover, some of its one-quark 
excitations are gapless quasiparticles, and the quark number density in the ground state is not zero. (In comparision, 
in the early investigations of the GN model at /z/ = and yU ^ the chirally broken phase was also found, but with 
zero quark number density and gapped quasiparticles [l^.) The region corresponding to this phase in the (^/,/i)- 
plane is not restricted. In contrast, in the chirally symmetric NJL model, i.e. with zero current quark mass, the 
(/i/, /i)-phase diagram does not contain the normal quark matter phase with such properties (i.e. with nonzero quark 
number density, gapless quasiparticles, and broken chiral invariance), at least at /i/ > (isj . 

Thirdly, we have found a rather rich quasiparticle excitation spectrum for the chirally symmetric phase 1. In this 
phase, if /i > z^, all quasiparticles are gapless. However, ii fi < v, then the gap is absent in the spectrum of only 
M-quasiparticles. 

We hope that our investigation of the GN phase diagram will shed some new light on the phase structure of QCD 
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at nonzero baryonic and isotopic densities. Thus, even in the most simple approach to the GN phase diagram we 
have found a variety of phases with rather rich dynamical contents. Obviously, a more realistic imitation of the QCD 
phase diagram requires to include also a nonzero bare quark mass, i.e. to study massive GN models, as well as to 
take into account the possibility of inhomogeneous condensates. 
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